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Abstract 

Two new one-parameter tracking behavior dark energy representations lo = a>o/(l + z) and 
uj = ujQe z ^ l+zS} I '(1 + z) are used to probe the geometry of the Universe and the property of dark 
energy. The combined type la supernova (SN la), Sloan Digital Sky Survey (SDSS), and Wilkinson 
Microwave Anisotropy Probe (WMAP) data indicate that the Universe is almost spatially flat 
and that dark energy contributes about 72% of the matter content of the present universe. The 
observational data also tell us that w(0) ~ — 1. It is argued that the current observational data can 
hardly distinguish different dark energy models to the zeroth order. The transition redshift when 
the expansion of the Universe changed from deceleration phase to acceleration phase is around 
zt ~ 0.6 by using our one-parameter dark energy models. 
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I. INTRODUCTION 

The SN la data suggest that the Universe is dominated by dark energy . Since 

1998, many dark energy models have been proposed in the literature. The simplest dark 
energy model is the cosmological constant model. However, the smallness of the value of 
the observed cosmological constant has puzzled theoretical physicists for a long time. For a 
review of dark energy models, see, for example references |^| and jsj. Although there exist 
a lot of dark energy models, we are still not able to decide which model gives us the right 
answer and find out the nature of dark energy. From theoretical point of view, perhaps the 
lack of understanding of quantum gravity is the main reason. To advance our understanding 
of dark energy, we may use observational data to probe the nature of dark energy. It is not 
practical to test every single dark energy model by using the observational data. Therefore, 
a model independent probe of dark energy is one of the best choices to study the nature of 
dark energy. 

The usual model independent method is through parameterizing dark energy or the equa- 
tion of state parameter u(z) of dark energy. The simplest method is parameterizing w(z) 



constant. To model the dynamical evolution o 



as the power law expansion u>(z) = 2 



N 



i=0 



LO;Z" 



dark energy, we can parameterize uj(z) 



. -Recently a simple two-parameter 



11 



model u){z) = ujq + u> a z/(l + z) was extensively discussed [1C , 
Padmanabhan later modified this two-parameter model as u(z) = ujq + 
More complicated forms of uj(z) were also discussed in the literature 



12 . Il^ . Jassal, Bagla and 




191. 



Instead of parameterizing oj(z), we can also parameterize the dark en ergy density itself, like 



a simple power law expansion fl 

26 



constant parameterization 



27 



21 



22 



28 



23, y, 



23J,|2J,|2J 

and the piecewise 



In [15[ , Gong used the SN la data to discuss some 
two-parameter representations of dark energy in a spatially flat cosmology. It was found that 
the SN la data marginally favored a phantom-like dark energy model. It was also found that 
the transition redshift zt ~ 0.3. In this paper, we propose two one-parameter dark energy 
models uo = uq/(1 + z) and uo = uo e z ^ 1+z ^ / (1 + z) and we use the SN la, the SDSS and the 
WMAP data to probe the geometry of the Universe. We also compare these two models 
with two two-parameter dark energy models. 

This paper is organized as follows. In section II, We first use the ACDM model as an 
example to show the method of fitting the whole 157 gold sample of SN la data compiled 
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in feci, the parameter A measured from the SDSS data |31| and the shift parameter 1Z 
measured from the WMAP data |26l . I27| to a dark energy model. The parameter A is a dark 
energy model independent parameter found by Eisenstein et al. in [3l| when they analyzed 
the large scale correlation function of a large spectroscopic sample of luminous, red galaxies 
for SDSS. It is related to the path from z = to z = 0.35. 1Z is the shift of the positions of 
the acoustic peaks in the angular power spectrum due to the effect of changing the values of 
Q m o and Qko on the Cosmic Microwave Background anisotropy. In section III, we propose 
two new tracking behavior one-parameter dark energy representations uo = loq/{\ + z) and 
u = ujQe z ^ l+z ' > I '(1 + z). In section IV, we fit the models to the observational data. In section 
V, we fit two two-parameter parameterizations u = u> +uj a z/(l+z) and u> = u Q +u a z/ (1+z) 2 
to the observational data. In section VI, we conclude the paper with some discussion. 



II. ACDM MODEL WITH CURVATURE 



For the simplest ACDM model where the dark energy is the cosmological constant, i.e., 
p = —p = A, we have 



h 2 = H 2 [n m0 (i + zf + n r0 (i + z) 4 - n k0 (i + zf + 1 + n k0 -n m0 -n 



rO 



where Q m (Q r ) = 8irGp m (p r )/3i?Q, = 8-35 x 10 5 |32|, and Q k = kja 2 H'^. The parame- 
ters fl m o and Qko are determined by minimizing 

[/WO - M*)] 2 , (A-0A69) 2 , (ft-1.716) 2 



x 2 = E 



(2) 



af 0.017 2 0.062 2 

where the extinction-corrected distance modulus f/,(z) = 51og 10 (dL(- 2; )/Mpc) + 25, the lumi 
nosity distance is 

'life 



di,(z) = a (l + z)r(z) = — — ^— Z ' sum 



dz' 



a H Jo E(z') 
dz' 



'\k\ 

-—7== sinn \ \/\ n ko\ , ... ,, 

H ox ftn^ r J ° E{z ^ 



(3) 



smn( 

E(z) 



k\x)/yj\k\ = sin(x), x, sinh(x) if k = 1, 0, — 1, the dimensionless Hubble parameter 
H(z)/H , the parameter A is defined as |3J| 



A 



rnO 



0.35 



0.35 



^(0.35) \n k0 



■sinn 



\Q 



k0\ 



0-35 dz 



E(z). 



1/3 



0.469 ±0.017, 



(4) 



3 



a 



0.1 



0.08 



0.06 



0.04 



0.02 



-0.02 



-0.04 



-0.06 



-0.08 



0.2 




0.22 



0.24 



0.26 



0.28 



0.3 



0.32 



0.34 



0.36 



0.38 



mO 



FIG. 1: The la, 2a and 3a contour plots of f2 m o and Qko for the model with the cosmological 
constant as dark energy 



the shift parameter 
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27] 



K 



m0 



i Q 



zSinn 



n 



n 



feO I 



k() I 



z is dz 



o E(z) 



1.716 ±0.062, 



(5) 



zi s = 1089 ± 1 [32J and Oi is the total uncertainty in the SN la data. In other words, we use 
the 1 57 g old sample SN la data compiled in jsij, the parameter A measured from the SDSS 
data [3III and the shift parameter 1Z measured from the WMAP data (26, 27] to find out the 
parameters Q m o and Qko- The nuisance parameter Hq appeared in Eq. (jSJ) is marginalized 
over with a flat prior assumption. Since Hq appears linearly in the form of 5 log 10 Hq in 
X 2 , so the marginalization by integrating £ = exp(— x 2 /2) over all possible values of H 
is equivalent to finding the value of Ho which minimizes \ 2 if we a l so include the suitable 
integration constant and measure function. 

The best fit parameters to the combined SN la, SDSS and WMAP data are f2 m0 = 
0.28 ± 0.03 and Q k0 = 0.004 ± 0.04 with X 2 = 177.14. The contour plot of Q m0 and Q k0 is 
shown in Fig. 
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III. ONE-PARAMETER PARAMETERIZATION 



The main goal of this work is to study the geometry of the Universe and the property 
of dark energy by using observational data. In the Introduction, we mentioned several 
different parameterizations. Those parameterizations have two or more parameters, so it is 
difficult to use those parameterizations to investigate the geometry of the Universe because 
we have to add two more cosmological parameters Q m o and Q^o to the model. Therefore, 
it is better that the parameterization has only one parameter so that the whole model has 
three cosmological parameters. Of course, this kind of parameterization assumed that the 
dark energy has evolution and it is not a cosmological constant. 

Model 1: To make the parameterization physical, we look for tracking behavior represen- 
tation, i.e., we require uj(z = oo) = 0. We first consider a simple one-parameter dark energy 
representation 

«M = ^ (6, 

At early times, z ^> 1, lu(z) ~ 0. In the far future, 1 + z — > and U)(z) — ► — oo. This 
simple parameterization has future singularity. The energy conservation equation of the 
dark energy can be written as 

dlnp = 3(l + a;) . . 

dz 1 + z ' U 

The acceleration equation is 

a 4ttG 



a 3 
Combining Eqs. (jOJ and (J2J), we get 



(p m + 2p r + p + 3p). (8) 



p = p (l + ,) 3 exp(^). (9) 

It is obvious that p ~ e 3wo po(l + z) 3 when z ^> 1 and p — > oo when z — > — 1. At early times, 
the energy density looks like matter with effective Q m o = e 3uJ °Qo. Here Q = 8ttGp/3Hq and 
Qo = 1 + Qko — H m o — Q r Q. This model may be thought as a unified model of dark matter 
and dark energy. The sound speed of the model is 

2= dp 2uj (1 + z)+3uj 2 

s dp 3(l + zy + 3ooo(l + z)- 1 ' 

So c 2 s0 = (2oo + 3wj})/(3 + 3lo ) ^ 0. When 1 + z + oo > and 1 + 2 + 1.5cj < 0, c^ > 0. 
For any unified theory of dark energy and dark matter, it was shown that the wave number 



k dependence of density perturbation growth due to the presence of a nonzero sound speed 
for a period of time produces unphysical oscillations or exponential blow-up in the matter 
power spectrum |33l . |3J, |35(. Therefore this model as the unified model of dark matter and 
dark energy is not feasible. The model with interactions between dark energy and dark 
matter was discussed in j^fj. 

Now we consider the model as a dark energy model. As we saw above, the dark energy 
behaves as ordinary matter at early times, we can interpret this as the tracking behavior, 
i.e., the dark energy tracked the matter at early times. The total effective matter density 
is ttf l0 = tt m0 + e 3a; °fi , and we expect that e 3w °fi < H m0 . Substitute Eqs. ® and © 
into Eq. (jSJ) and neglect the radiation contribution, we find that the transition redshift z T 
satisfies the following equation 



:i + ^) 3 



0. (11) 



Model 2: Now we consider another one parameter dark energy parameterization 

u{z) = J^_ e */d+*). (12) 
v ' 1 + z 

For this model, uj(z) ~ when z> 1. The major difference between this model (fT2*|) and 
the model © is that w(z) — ► as z — > —1 for this model. These two models have almost 
the same behavior in the past and very different behavior in the future. Combining Eqs. 
(J7j) and 0, we get 

p = p (l + zf exp (3cu e 2/(1+2) - 3cu ) . (13) 

It is obvious that p ~ e 3uj ° e ~ 3uJ ° p (l + z) 3 when z>l and p ~ e~ 3uJO p (l + z) 3 when z — > — 1. 
At early times, the energy density looks like matter with effective Q m o = e 3aJ ° e ~ 3a;o no, and 
it behaves like matter with effective Q m Q = e~ 3uJ °Qo in the far future too. This model may 
also be thought as a unified model of dark matter and dark energy. The sound speed of the 
model is 

2 _ dp _ 2cj (1 + z)e z ^ l+ ^ + cu (l + 3cj e 2 /( 1+2 ))e 2 /( 1+2 ) 
° s ~ dp~ ~ 3(1 + zf + 3c^/( 1+2 ) ' ( ' 

So c 2 s0 = ujq < and it also produces exponential blow-up in the matter power spectrum. 
This model as the unified model of dark matter and dark energy is not feasible too. Again 
we consider this model as a dark energy model. One key feature of this model is that the 
model behaves like matter both in the past and the future. The Universe will expand with 
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deceleration in the future. In the past, the dark energy tracked the matter. The total 
effective matter density is Q^ = Q m o + e 3uJ0 ^ e ~^Q , and we expect that e 3uJo( - e ~^Q "C fi m o- 
Substitute Eqs. ()12|1 and (jT3j) into Eq. (jSJ) and neglect the radiation contribution, we find 
that the transition redshift zt satisfies the following equation 

, z 



1+z) 



n 



mO 



3^Q 

I + — — exp( 



1 + r 



fioexp (3cj (e 2/(1+2) - 1 



0. 



(15) 



IV. DATA FITTING RESULTS 



By fitting the model 1 to the combined SN la, SDSS and WMAP data, we get f2 m o = 
0.25 ±0.05, Q k0 = -0.009 ±0.05 and u = -1.1 ±0.2 with X 2 = 175.4. If we take the model 
1 as a unified model of dark energy and dark matter, we find that the best fit parameters to 
the combined SN la, SDSS and WMAP data are tt k0 = -0.05 ±0.04 and u = -0.42 ±0.04 
with x 2 — 203.6. Since Ax 2 = 203.6 — 175.4 = 28.2, we conclude that this model as a unified 
model of dark matter and dark energy is not a viable model. The contour plots of fl m o and 
Qko by fixing ujq at its best fit value —1.1 are shown in Fig. El The contour plots of Q m o 
and ujq by fixing Q k0 at its best fit value —0.009 are shown in Fig. El The evolution of uj(z) 
is shown in Fig. |H1 Substitute the best fit values to Eq. (JTTJ), we get the transition redshift 

= 0.56. The results are summarized in Table HJ 

If we use SN la only, then we get fl m0 = 0. 331^12, ^fco — ~0.33lil'll and co = — 7.5i 3 3 8 
with x 2 = 171.9. 

By fitting the model 2 to the combined SN la, SDSS and WMAP data, we get fl m o = 
0.28 ±0.04, Q k0 = -0.0011^ and u = -0.97±° H with X 2 = 176.5. If we take the model 2 
as a unified model of dark energy and dark matter, the best fit parameters to the combined 
SN la, SDSS and WMAP data are tt k0 = -0.10 ±0.05 and cu = -0.22±g;g§ with X 2 = 233.2. 
Again this model as a unified model of dark matter and dark energy can be firmly ruled out. 
The contour plots of Q m Q and Q k o by fixing ujq at its best fit value —0.97 are shown in Fig. 
0] The contour plots of Q m Q and ujq by fixing Q k0 at its best fit value —0.001 are shown in 
Fig. El The evolution of u(z) is shown in Fig. El Substitute the best fit values to Eq. (JEJ), 
we get the transition redshift zt = 0.66. The results are also shown in Table HI 

If we use SN la only, then we get O m0 = 0.33_o.20) ^fco = — 0. 341^32 and too = — 7.4^3 3 8 4 
with x 2 = 171.9. 

Although the two models discussed in the previous section have very different future 
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FIG. 2: The la, 2a and 3a contour plots of f2 m o and Qko for the parameterization u = ojq/(1 + z). 
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FIG. 3: The la, 2a and 3a contour plots of f2 m o and ojq for the parameterization uj = loq/{1 + z). 
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FIG. 4: The la, 2a and 3a contour plots of f2 m o and fifco for the parameterization cu = u>o exp(z/(l+ 
z))/(l + z). 
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FIG. 5: The la, 2a and 3a contour plots of Vt m Q and Uq for the parameterization oj = uq exp(z/(l + 
z))/(l + z). 



behavior, they both fit the current data as well as the ACDM model. This may suggest 
that the current data fitting method cannot distinguish models with very different future 
behavior. The best fit results also show that the Universe is almost spatially flat, and that 
the best fit results using the combined SN la, SDSS and WMAP data are different from 
those using SN la alone. 



V. TWO-PARAMETER PARAMETERIZATION 



In this section, we consider spatially flat cosmology onh, 
Model 3: We first consider the parameterization 



gy only. 



uj = uj + — ^— . (16) 

When z 1, we have to ~ ujq + ui a . u — > ±oo when z —* —1. Combining Eqs. (J7J) and (fTfij). 
we get the dark energy density 

n = Q (l + zfi^+^a) eX p(-3cj a z/(l + z)), (17) 

where Qo = 1 — O m0 — Q r o. Substitute Eqs. (|T5|) and ([T7jl into Eq. (jHJ) and neglect the 
radiation contribution, we find that satisfies the following equation 



^mO + (1 ~~ ^mO, 



(l + 3^o + (1 + exp (Z^) = 0. (18) 



h0.35 



The best fit to the combined SN la, SDSS and WMAP data gives that u = -1.13Io;b 6) 
uj a = 0.95l?;95 and VL m0 = 0.28 ± 0.04 with x 2 = 175.62. Substitute the best fit parameters 
into Eq. (|18|). we get zt = 0.56. The results are summarized in Table The contour plots 
of lvq and u a by fixing Q m0 at its best fit value 0.28 are shown in Fig. |U1 The evolution of 
u)(z) is shown in Fig. IHl 

n 

Model 4: Next we consider the following parameterization [14], 

u = UJo + JiTW- (19) 

When z ^> 1, we have to ~ loq. When z — > — 1, we have u; — > ±oo. Substitute Eq. (fTTJj) into 
Eq. (JZJ), we get the dark energy density 

n(z) = fi (l + z) 3(1+wo) exp(3cj a 2 2 /2(l + ^) 2 )- (20) 
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FIG. 6: The contour plot of ujq and u a for the parameterization uj = ujq + u a z/(l + z) 

Substitute the above two equations (|T5|) and (J2Uj) into Eq. (JHJ) and neglect the radiation 
contribution, we find that z T satisfies the following equation 

n m0 + (1 - Q m0 ) (l + 3u + 7f^) (i + exp i^^f) = °- ( 21 ) 

The best fit to the combined SN la, SDSS and WMAP data gives u = —1.51 ± 0.6, 
ijj a = an d ^mo = 0.27 ±0.04 with x 2 = 174.0. Substitute the best fit parameters into 

Eq. (121}) . we get z^ = 0.38. The results are summarized in Table |1J The contour plots of cuo 
and io a by fixing Q m0 at its best fit value 0.25 are shown in Fig. [7| The evolution of uj(z) is 
shown in Fig. |S1 

VI. DISCUSSION 

We discussed two one-parameter dark energy parameterizations and two two-parameter 
dark energy parameterizations. For the two one-parameter dark energy parameterizations, 
we consider curved cosmology, so that we have a total of three cosmological parameters: 
^mo> ^fco an d ^o- F° r the two two-parameter dark energy parameterizations, we consider 
flat cosmology only, again there are three cosmological parameters: fl m o, an d ui a . These 
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FIG. 7: The contour plot of ujq and uo a for the parameterization to = luq + u a z/(l + z) 2 



TABLE I: Summary of the best fit parameters. 



Model 






UJ 






x 2 


1 


0.25 ±0.05 


-0.009 ±0.050 


-1.1 ±0.2 


N/A 


0.56 


175.4 


2 


0.28 ±0.04 


-n noi +0 - 046 

U.UUl_o.0 45 


u - y/ -0.19 


N/A 


0.66 


176.5 


3 


0.28 ±0.04 


N/A 


1 1O+0.35 
L - LO -0.26 


u - yo -1.95 


0.56 


175.62 


4 


0.27 ±0.04 


N/A 


-1.51 ±0.6 


A q+3.8 
4 "J-4.8 


0.38 


174.0 



different classes of three cosmological parameter models fit the observational data almost 
equally well because they have almost the same minimum value of \ 2 as shown in table UJ 
However, they have very different behavior. At early times, the Universe is dominated by 
matter or radiation, the dark energy is subdominant, so the contribution of dark energy 
to the background evolution is not important and the data may not be used to distinguish 
the early behavior of dark energy to the 0-th order. From Fig. |H1 we see that the future 
behavior of uj(z) are also very different. For the model 1, lu(z) — > — oo in the future. For 
the model 2, to{z) ~ in the future. For the model 3, uj(z) — > — oo in the future. For 
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Redshift z 



FIG. 8: The behavior of lo{z). The solid lines plot oj{z) by using the best fit parameters and 
the dash dotted lines are for la errors. For the two-parameter representations, the la errors are 
obtained from the la contours of ujq and tu a . 

the model 4, lu(z) — > +00 in the future. So the data may not used to distinguish the 
future behavior of dark energy to the 0-th order either. We need to invoke at least linear 



perturbation method to discuss dark energy models. In [37[ , the authors use the concept of 
the minimal anti-trapped surface or the assumption that the energy momentum content of 
the observable Universe does not change significantly in comoving coordinates to study the 
fate of our universe. They found that it is impossible to confirm the accelerating expansion 
with current observed dark energy value Qq ~ 0.7 if the dark energy is not a phantom. These 
results more or less support our conclusion. The dark energy in the models 1 and 2 tracked 
the matter in the past. The two models both suggest that the Universe is almost spatially 
flat. All the models suggest that zt ~ 0.6 and u)(0) ~ —1. These results are consistent with 
those derived from the simplest ACDM model. However, it was found that z-r ~ 0.3 by using 
SN la data only [l^. So the results by using combined SN la, SDSS and WMAP data are 
different from those by using SN la data only. More thorough studies are needed to make 
more concise conclusion. Finally, we would like to mention that the one and two-parameter 
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representations of dark energy models have their own limitations as discussed in |38j . These 
parameterizations do not accommodate the possibility of rapid evolution of dark energy. 
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